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1. Introduction 

Entanglement is a property which is central to quantum mechanics, with bipartite 
entanglement being readily producible experimentally. There has been some progress 
in the production of tripartite entangled beams, although the entanglement is often 
obtained by mixing squeezed vacua with linear optical elements PlEj- Other methods 
which create the entanglement in the actual nonlinear interaction have been proposed, 
using both cascaded and concurrent x*"^"* processes jSl IH Ej- In this article we 
investigate the fundamental limits to the achievable tripartite entanglement available 
from two processes which utilise concurrent and cascaded nonlinearities. We evaluate 
the continuous variable tripartite entanglement criteria of van Loock and Furusawa jB] 
for two different interaction Hamiltonians and then calculate the performance of the 
corresponding intracavity systems which contain the interactions described by these 
Hamiltonians. 

2. Criteria for tripartite entanglement 

We will first describe the inequalities which should be violated to demonstrate that a 
system demonstrates true continuous variable tripartite entanglement. For three modes 
described by the annihilation operators dj, where j = 1,2,3, we define quadrature 
operators for each mode as 

Xj = dj + aj, Yj = -i{dj - a]), (1) 

so that the Heisenberg uncertainty principle requires V{Xj)V{Yj) > 1. A set of 
conditions which are sufficient to demonstrate tripartite entanglement have been derived 
by van Loock and Furusawa [Hj, without making any assumptions about Gaussian 
statistics. This is in contrast to the usual conditions used to demonstrate continuous 
variable bipartite entanglement, which were developed by Duan et al using the properties 
of the covariance matrix for Gaussian variables [7]. Using our quadrature definitions, 
these conditions give a set of inequalities, 

^2 = V{X^ - X2) + V{Y^ + % + %) >4, 

\/i3 = V{X^ - X3) + V{Y, + % + %) >4, (2) 

^23 = V{X2 - X3) + V{Yi +% + %) >4, 

where ^(^4) = {A^) — {A)"^. If any two of these inequalities are violated, tripartite 
entanglement is guaranteed. 

We note here that, although states which violate these inequalities are sometimes 
called continuous variable GHZ states, actual GHZ states rely on perfect correlations 
and are not at all of a statistical nature jHl E] • This terminology would be accurate if 
the correlations above actually vanished, so that the states involved were eigenstates 
of amplitude quadrature differences and phase quadrature sums. This would, however, 
require infinite squeezing of the electromagnetic field which in turn would require infinite 
energy. Although this situation has been predicted in some works as the result of 
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linearised fluctuation analyses, it always occurs in regions where this form of analysis is 
not valid. Hence we will not use this terminology in this article. 

3. Travelling wave models 

In this section we will examine simplified models of two processes which utilise two 
concurrent or cascaded nonlinearities to produce tripartite entanglement. We will 
proceed both by solving stochastic equations which are equivalent to the Heisenberg 
equations of motion derived from the interaction Hamiltonians and also approximate 
operator equations obtained by ignoring pump depletion. While these procedures are not 
intended to model realistic physical processes with propagating waves, which would need 
a more complicated treatment PU] ? they do give some idea of the degree of entanglement 
available from these Hamiltonians. In the next section we will perform more realistic 
analyses of the same processes inside pumped optical cavities. 

3.1. Two cascaded nonlinearities 

Ferraro et al have proposed a linked nonlinear process which links five different modes jl] 
and was first investigated in terms of photon statistics by Smithers and Lu It has 
recently been analysed in terms of its nonlocal properties by Ferraro and Paris ^2] , with 
the theoretical analyses of all these works using the undepleted pump approximation. 
However, it is known that the squeezing and bipartite entanglement in travelling-wave 
T^(2) processes do not become perfect as the interaction strength increases, but reach 
some finite limit and then decrease ^21 ^1 El- We will therefore quantise the pump 
modes and analyse this system using the full stochastic equations of motion, which 
must be done numerically. The process involves modes which we will describe using the 
operators fei(a;i), 62('^2), ai('^3), ^^2(1-^4) and a3(ci;5) where the frequencies obey 

The interaction Hamiltonian can then be written 



which we see describes the same process as ref. |3] once all the interacting fields are 
quantised. Note that we have changed the indices of modes 2 and 3 compared with 
those used by Ferraro et al so that both our 1 and 2 are produced by the xi interaction. 
This Hamiltonian describes a downconversion process cascaded with a sum-frequency 
generation process where one of the downconverted modes becomes a pump mode for 
the frequency generation process. 

Before we develop and investigate the full equations from (jH), it is instructive to 
examine the analytical solutions which may be obtained using an undepleted pumps 
approximation. Setting k,i = Xi(^i(0)) and K2 = X2(^2(0)) as real positive constants, we 




(4) 



^202. 

dt 
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find the Heisenberg equations of motion, 

dai ^1 
— = .,a2, 

^ = Kia\ - K,2a3, (5) 
dt 

dOg 

dt 

We find that there are two classes of solutions, depending on whether > or i^l < nl. 
In the first case, where k| > k^, we set = ^y — to find 

„ ,, Ko — K?cosr2t„ , , Kisinfit^t, , niHoicosQt — 1) , , 

= ^—^2 + ^ 4(0) + ^ ^2 ^4(0), 

„ , , Kisinfit^t, , ^ , . ^ K2sinr2t^ , , 

a2{t) = ^-^ a{{0) + d2{0) cos fit - ^-^ 03(0), (6) 

„ KiK2(l — cos fit) ^+ , , K2Sinfit„ , , KoCOsfit — , , 

«3W = ^2 ^«l(0) + ^^^2(0) + -^^^ ^63(0), 

which, beginning with all the a modes initially as vacuum, gives the solutions for the 
intensities as 

i sin^ fit 
(aW) = ^ ) 

. _ nl4{cosQt-lf (7) 
<;a3«3/ ^ , 

{a{ai) = (0^02) + (403)- 
We see that these are the same as the analytical solutions given by Ferraro et al jl], 
once the change of indices is taken into account. 

We can use the solutions of © to find expressions for the quadrature variances and 
covariances, 

4k^k|(1 - cos fit) - 2Kf sin^ fit 



V{X^) = y(Fi) = 1 + 



fi4 

2 „;„2 



^ (8) 

V^(Xi,X2) = -ViY,,Y2) = ^^1^ (4 - «:^osfit) , 

V{Xi,X3) = -V{Yi, Y3) = [f^l cos' nt + Kl- {kI + 4) COS fit] , 

ViX2, X3) = ViY2, Y,) = sin fit (1 - cos fit) . 

These expressions contain all the information necessary to express the van Loock- 
Furusawa correlations, since V{Xi-Xj) = V{Xi) + V{Xj)-2V{Xi, Xj) and V{Yi + Y2 + 
Y3) = V{Y,) + + + 2{V{Y,, %) + %) + %)). The expressions 

for the van Loock-Furusawa correlations obtained by combining these variances and 
covariances are shown in figure ^ for K2 = 1.8ki. We see that this Hamiltonian provides 
tripartite entanglement over a range of scaled interaction time. 




Figure 1. The analytical solutions of the van Loock-Furusawa correlations for the 
Ferarro scheme, with K2 — 1.8ki. The line at 4 represents the level beneath which at 
least two of the correlations must fall to exhibit tripartite entanglement. 



If > K2 the solutions are not periodic so that the undepleted pumps 
approximation is of hmited vahdity, being expected to give accurate answers only for 
short times. However, setting = — we find 



K|coshCt-KL , '^isinhCt.t/n\ , '«i«;2(l - coshCt) .t, 
^2 «i(0) + ^ «2(0) + ^ Oaf 



, , Ki sinh t / N / N , > K,2 sinh Ct , , 
Ht) = ^ alio) + 02(0) cosh a - ^ a3(0), (9) 

. KiK2{coshCt-l) KasinhCt. , K^-K^coshCt. 
flsW = ^2 "1(0) + ^ "2(0) + 03(0), 

with the mean intensities being, again with these modes beginning as vacuum. 



Kfsinh 



K?K^(coshCt - 1)^ (10) 



^4 

{a[ai) = (4^2) + (4^3)- 

We note here that these solutions have previously been given by Smithers and Lu JT] 
and that we have presented them here because we intend to use them to find analytical 
expressions for the correlation functions of interest. In this case the solutions for the 
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Figure 2. The analytical solutions of the van Loock-Furusawa correlations for the 
Ferarro scheme, with ni — \.1ki. The line at 4 represents the level beneath which at 
least two of the correlations must fall to exhibit tripartite entanglement. 



variances and covariances are found as 

2k\ \k\ sinh^ C^t^ k\{2-2 cosh (t)] 



V{X^) = = 1 

V{X2) = V{Y2) = 1 
V{Xs) = V{Y,) = 1 



^4 5 

2kI sinh^ 

2fi;^K|(coshCt - 1)2 

(11) 



^4 



V{X^,Xs) = -V(Y„Y,) = ^ [{^ + - cosh a) + /c^cosh^a] 



r(X2, Xs) = %) = — ^ sinha (cosh a - 1) • 



The system also exhibits tripartite entanglement in this regime, as can be seen in figureEl 
However, as the solutions for the intensities are hyperbolic, they quickly increase to the 
point where the undepleted pumps approximation will lose its validity. Neither of the 
analytic treatments used above is useful for the case where for which we will 

employ stochastic integration. 

In developing our full equations of motion, we will follow the approach of Huttner 
et al. ^ni; (see also Caves and Crouch ^7]) treating the interacting fields in terms of 
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the photon fluxes rather than in terms of energy densities. As stated in Ref. ^Hl, this 
approach avoids problems which could arise, especially with the quantisation volume, 
if we were to work with the normal Hamiltonian approach. With the appropriate 
momentum-space operators, we use the well-known mapping onto stochastic differential 
equations in the positive-P representation ^H] to calculate the development of the fields 
as they traverse the medium. We note here that this phase-space representation allows 
for an exact and complete mapping of our Hamiltonian onto stochastic differential 
equations. We consider here the case of one dimensional propagation, which is valid 
for the case of colinear pumping. We also note here that this approach assumes that 
the medium is not dispersive for the interacting fields, which is a difficult condition to 
meet with existing materials. In this approach, the operator 

N{zo, UJm) = a^(2o, UJm)a{Zo, UJm), (12) 

for example, is the number operator for photons at frequency Um which pass through 
a plane at z = zq during a chosen time interval. The operators a'^{z,ujm) and a{z,Um) 
then obey bosonic spatial commutation relations, 

[a{z,uJi),a\z',ujj)] = 6ij6{z - z'), (13) 

and similarly for the bj operators. The nonlinear momentum operator for this system 
is found as 

Gni{z) = ih (^Xib\aia2 + X2blald3^ + h.c. (14) 

As shown by Shen ^H], we can write an equation of motion for the density matrix of 
the system, 

i^^= [p{z),GU^)] , (15) 

which allows for the calculation of steady-state propagation, exactly as required for 
continuous pumping. Physically, the density matrix, p{z), describes an ensemble of 
steady-state systems which has all the statistical properties of the fields at point z. 
Equation (fT3j) provides a full description of the interacting fields of our model, but is 
extremely difficult to solve directly. 



Continuous variable tripartite entanglement from twin nonlinearities 



8 



Therefore, following the standard procedures we map the master equation onto 
the Fokker-Planck equation for the positive-P pseudoprobability distribution, 



dP 

d2 



d 

dai 



d 



dai 



d 



d 



da2 

d 
das 

d 



d 

X20i2(^2 + -7r^X2atl32 



dai 



d 



dP2 



-Xl"l"2j 



{-X2(y2(yt) 



d 



d 



{-Xi^tat) 
- {-X2ata3) 



(16) 



1 

+ 2 



daida2 



2xiPi 



-2X2(^3 - 



q2 

da^da2'^^^^^ 
:2X2a3 



Pia.z) 



da2dP2 da2dP2 

where a = {ai,af ,a2,a2 ,a3, a^ , Pi, Pi , (32, Pt)- always with the positive-P 

representation, stochastic averages of products of the variables represent normally 
ordered operator expectation values, with there being correspondences between 
j',Pj,Pj' and aj,d], bj,b]. We now map this Fokker-Planck equation onto the 
following set of stochastic differential equations in Ito calculus, 

dai 



dz 
daf 
~d^ 
da2 
~d^ 

da^ 

~di 

das 

dz 
dc4 

dz 
dPi 

dz " 
dPl 

dz 
dP2 

dz " 

dPl 
dz 



X 



Xia2pt + \l iV2 + m) 



XiafPi - X2asPt 

- Xi(^il3t - X2atP2 

X20i2p2, 

- X2cytp2^ 

-Xi"ia2, 
-X20i2at + 



XiP 



XiPt 



-{fli 



X20i3 



{V2 - m) - 



X2«3 



(^8 - iVe) 



[17) 



X20:3 



iV7 + iVb] 



-X2a2 «3 + 



X2(y3 
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5.5 




Figure 3. The scheme of reference ^ with everything symmetric, /3i(0) ~ /?2(0) — 10^ 
and Xi = X2 = 10~^, averaged over 1.06 x 10^ stochastic trajectories. The horizontal 
axis is the scaled interaction length, ^ = |/3o|x-2^- The line at 4 represents the upper 
bound for true tripartite entanglement. 

which we may solve using stochastic integration. The real Gaussian noise terms have 
the correlations 

r]j{z) = 0, r]j{z)r]k{z') = 6jkS{z - z'). (18) 

The results of stochastic integration in two different parameter regimes are 
presented in figure El figure |3] and figure El All show genuine tripartite entanglement 
over some interaction range and all begin with the output modes as vacuum. For the 
parameters used in figure El with XiPii^) = X2/32(0), the output intensities increase 
monotonically over the range shown and the entanglement disappears. This is unlike 
the situation of figure (H where X2/32(0) = 2xi/3i(0), and both the output fields and the 
van-Loock Furusawa correlations oscillate over a short interaction length. Although the 
field intensities for the first situation obviously cannot increase indefinitely, this contrast 
between monotonically increasing and periodic behaviour of the intensities, depending 
on the ratios of the pumping and interaction strengths, was mentioned by Smithers and 
Lu [HI- 

3.2. Two concurrent nonlinearities 

Another possibility for a travelling wave model is to have a single crystal with two 
concurrent nonlinearities, each pumped by different modes. This could be achieved 
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Figure 5. The intensities produced by the interaction of the scheme of reference 0] 
for the same parameters as used in figure 01 
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either with different polarisations or different frequencies. In this section we will 
consider a crystal which is pumped at frequencies Ui and 1^2 to produce modes at 
uj3,uj4: and UJ5, where ui = + uJi^ and (^2 = 1^4 + 1^5- With all modes quantised with 
the operators 61 (a;i), 62(^^2), 0^1(^^3)5 0^2(^4), ^(i-^s), the interaction Hamiltonian for this 
scheme becomes 

-f^int = i^ (^Xi^i«i4 + X2^2a24) + h-c, (19) 
where the Xj represent the nonlinear interactions. 

We will first examine a simplified analytical model for the propagation of fields 
described by this Hamiltonian. Assuming perfect phase-matching and ignoring pump 
depletion, we may set 71 = Xi(^i(0)) and 72 = X2(^2(0)), where the pump fields are 
initially intense coherent states, to find the Heisenberg equations of motion, 

dai . 

^ = lia\ + 72aL (20) 



Setting f2 = v^7i + 7I, we may solve these linear operator equations to find the solutions 

. /,x . ^2 + 7? cosh fit 7isinhfit . 7172 (cosh fit - 1) 
ai(t) = ai(0) — + a^(0) + 03(0) — , 

^ / N ^t/ .7isinhfit „ , ^ „t /„x 72 sinh fit 

a2(t) = d{{0)^^-^ + 02(0) coshfit + G^(0)^^-^^ , (21) 

- r+\ - /n^ ^1^2 (cosh fit - 1) 72sinhfit . 7f + 7|coshfit 
a3(t) = ai(0) — + a^(0) + 03(0) — , 

which contain all the information required to calculate the desired correlations. With 
all the output modes initially vacuum, we find the average intensities 

. 7^sinh^fit 
{a[ai) = — , 

(0^02) = sinh^ fit, (22) 

, 7|sinh^fit 
(aaQa) = ^ , 

with the variances and covariances for the quadratures being found as 

7isinh^fit 



fi2 

,2 , 



V{X2) = ¥{¥2) = 1 + 2 sinh^ fit, 

V(X,) = = 1 + 22i^. 

71 sinh 2fit (23) 



\/(Xi,X2) = -F(Fi,F2) 
V{X2,X,) = -V{Y2,Y3) 



fi 

27172 sinh^ fit 

fi2 ' 

72 sinh 2fit 
" fi ■ 
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We can now calculate the van Loock-Furusawa correlations, finding 




(24) 



which, for 71 = 72 = 7, simplify to give 

V12 = 1^23 = 5 + 9 sinh^ nt - 3^2 sinh 2nt, 

2 r (25) 

= 5 + 6 sinh^ nt - 2V2 sinh 2fit, 

from which we can see that tripartite entanglement is present for shortish interaction 
times, but does not increase with Qt. In fact, the minimum of V12 and V23 is found at 
ilt = cosh~^ \/2, and has a value of 2. 

We will now turn to stochastic integration of the full equations, without using the 
undepleted pump approximation. In what follows we will set the interaction strengths 
equal, Xi = X2 = Xi these are the conditions which give the maximum of tripartite 
entanglement for this sytem. We define the nonlinear momentum operator for this 
process as 



From this operator we may make a mapping onto a Fokker-Planck equation for 
the positive-P pseudoprobability distribution, from which we find the Ito stochastic 




(26) 
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differential equations, 



1^6 J 



dz 
da 



_ = xa+(32 + y — (^3 - 1^7) , (27) 
d/5i 



d^ 
d^ 

d^ 
dz 

d(3+ 



-Xttl«2, 

-X«2a3, 



dz 

where the real Gaussian noise terms have the correlations 



7]^{z) = 0, vA^)Vk{z') = 6,kS{z - z'). (28) 

As in the previous section, stochastic averages of products of the variables represent 
normally ordered operator expectation values, with there being correspondences between 
a j , , Pj , ^^d dj , d] , bj , 6] . 

The results of stochastic integration of ()27|) are shown in figure IHl for the tripartite 
entanglement criteria, for parameter values x = 10~^, Pi{0) = ^2(0) = 10^ and 
ai(0) = 02(0) = 03(0) = 0. We see that the correlations are not symmetric, but 
that tripartite entanglement is available. The output intensities are shown in figure [71 
from which we again see that the inequalities are violated for relatively weak fields. For 
this system, the analytic results which we presented above give the same results as the 
stochastic integration over the interaction range shown. 



4. Doubled intracavity nonlinearities 

Now that we have demonstrated that the two Hamiltonians of section IHl can produce 
entanglement, we will examine the more realistic physical situations where the processes 
described happen inside pumped optical cavities. As optical cavities can be tuned to 
be resonant with only a certain number of modes, a fuller description of the physics 
involved can be given in a simple manner. 
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Figure 6. Positive-P solution averaged over 2.68 x 10^ stocliastic trajectories for the 
tripartite entanglement criteria with the doubly concurrent Hamiltonian of section llOl 
The solid line is V12 and V23, while the dashed line is V13. The horizontal axis is the 
scaled interaction length, ^ ~ \Pq\xz. The line at 4 represents the upper bound for 
true tripartite entanglement. 




Figure 7. Output intensities for the same system and parameters as figure El 
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4.1. Cascaded nonlinearities in a cavity 

To consider the scheme described by the interaction Hamiltonian of (jH) inside a resonant 
pumped optical cavity, we must add terms to the Hamiltonian, so that 

H = i/int + -f^pump + -f^damp, (29) 

where i^int has the same form as the interaction Hamiltonian of (jH), and 
^^pump = (ei&l + £2^2) + h-c. 

i^damp = ^{12 ^'4 + E r^^l j + h.c, (30) 

where /c = 1,2,3 and j = 1,2. In the above, the are the classical pump amplitudes 
at the two input frequencies and the F^^, are bath operators. The field operators now 
refer to the intracavity fields. This now becomes equivalent to the scheme of Guo et 
al jHj, who analysed it using quantum Langevin equations in an undepleted pumps 
approximation. Following the same procedures used in deriving (fT7j) and making the 
usual zero-temperature Markovian bath approximation, we find the positive-P equations 
which give a full quantum description of this system, 

dai 



dt 



-7iai + Xitt2 A + 



- 7iai + Xi"2Pr + \\ -Tr~ (^^2 + 1114) 



dt ' ' V 2 

da2 
Id 



-72^2 + Xitti A - X2a3pi + 



dt 
dtts 



72^2 + - X2atP2 + \l {V2 - m) - \l -rr- iVs - iVa 

-73a3 + X2a2/?2, 

= -73"^ + X2at(3t, 
ei - Ki/3i - Xi"ia2, 



dt 
dot 

dt 
dA 

dt 

df3l 
dt 

— = 62- K2P2 - X2a2aJ + J (?77 + iV5) 



'^2l3t - X2aJ«3 + Y (^8 + ^Ve 

Note that the correlations of the noise terms are now in time rather than in the spatial 
variable. These equations can be integrated numerically in any parameter regime, 
including near to any critical points of the system. However, as we will not be concerned 
with the behaviour of this system in the neighbourhood of any critical points, we will 
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proceed via a linearised fluctuation analysis. This involves separating the variables of 
the positive-P equations into their mean-field steady-state solutions plus a fluctuating 



part, e.g. ai 



+ 6ai. Solving the classical equations of motion to find the steady- 



state solutions, we may then write equations for the fluctuations from which we can 
calculate the output spectral quantities of interest |2II. Neglecting the noise terms in 
^T\ . we find the following classical equations for the interacting fields, 

- -7iai -Fxi«2Pi) 
ei - Ki(3i - Xittia2, 

62 - K2P2 - X2«2«3' (32) 

where the 7^ (kj) are the cavity damping rates for the aj 

We see that one possible set of solutions to these equations is 

af = 0, (33) 

PT = ^^/S' (34) 
which are reminiscent of those found for the well-known OPO below threshold. To 
examine the stability of these solutions, we write the linearised equation for the 
fluctuations, 

dSa = ASadt + BdW, (35) 

where a = [6ai,6af ,6a2,Sa2 ,Sa3,6a^,6j3i,6j3i ,6/32,^132]^ ^ ^ is the matrix of the 
noise terms of |T7jl . but with the classical steady-state solutions used in place of the 
stochastic variables and dW is a vector of Wiener increments. The drift matrix is found 



dt 
da2 

'df 

da-i 
dt 

dt 
d(32 
dt 



as 



A = [A^ A2] 



(36) 



where 



Ai 



-71 








Xi/5f 











-7i 


xiiPir 
















-72 





-X2ip;r 













-72 





-X2/3I 








x2/3r 





-73 














X2{f3*2r 





-73 





















-xi{a*2r 





-xMir 














-X2ia*,r 








-X2a-2 











-X2a3 


-X2ia*2r 






(37) 
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and 



A, 

































-X2«3 







-X2iaiy' 



















X2{a*2r 






















-Ki 














-1^2 














-fi2 



(38) 



As long as none of the eigenvalues of this drift matrix have a positive real part, the 
solutions are stable and the linearised fluctuation analysis should be valid. Although 
general analytical expressions for the eigenvalues can be found, these are rather 
complicated. In the simplifying case below threshold where we set the output loss 
rates equal, 7j = 7 and kj = k, we find a degenerate eigenvalue which can be positive. 



A 



-7 



xl4 



K 



This eigenvalue sets a condition for the relative strengths of the pumps. 



2 2 
Xi^i 



22^ 22 

X2^2 <1 ■ 



When this condition is violated, the below threshold solutions are unstable. 

The spectral correlations are found in the normal manner, via the equation 



S{lj) = {A + iu)~^ D {A^ - luo)' 



(39) 



(40) 



(41) 



where D is the diffusion matrix with the steady-state values of the fields, and the 
standard input-output relationships [221 • We will denote the output spectral correlations 
equivalent to the Vij of (j2I) as Sij. Although we were able to find analytic expressions 
for the correlations, these were exremely complicated and gave little insight, therefore 
we have presented numerical results in figure |H1 Our numerical investigations over a 
range of parameters did not find any tripartite entanglement noticeably better than 
that presented. We found that when the pumping rates increase above the value of 
ej = "Jk/xj, the entanglement rapidly disappears. What is immediately visible is that 
the inequalities are not all violated equally, although true tripartite entanglement is 
demonstrated. We note that numerical investigations show that this system exhibits 
a range of behaviours, with self-pulsing type oscillations and possible bistability for 
particular parameter regimes, but that here we are only interested in its suitability as 
a source of tripartite entanglement and an investigation of these effects is outside the 
scope of this work. 
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4 6 
CO (units of y) 



10 



Figure 8. Tripartite entanglement criteria for the system of section ETI for k = j 

Xi=X2 = 10"^ and ei = £2 = 0.9jk/xi- 



4.2. Intra-cavity concurrent twin nonlinearities 

In a similar manner to the preceding section ()4.H1 . we will now investigate the 
performance of an intracavity version of the system described by (fT^ . We again find the 
classical equations by dropping the noise terms in the appropriate positive-P equations 
and solve these for the mean fields in the steady-state. The equations for the mean 
fields are 



da I 

'df 

da2 

'df 

da-i 

dt 
dA 

dt 
dfh 

dt 



-7a«3 + XOi*2p2, 

£2 - lbf32 - Xtt2a3, 



(42) 



where the 7a,b are the cavity damping rates for the appropriate modes and the ej are 
the classical pumping terms. Although we have chosen the simple case where both 
pump modes have the same loss rate, as do the three signal modes, this is not essential, 
although it does serve to simplify our analysis. We may now solve the above equations 
to find the steady-state mean field solutions and the conditions for the stability of the 
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Figure 9. Tripartite entanglement criteria for the system of section l?!^ for = Jb = 
1,X^ 10"^ and e = OJeth- 



linearised fluctuation analysis. 

To find analytical expressions, we will set €i = €2 = e. We find that there is an 
oscillation threshold at the value eth = 'Jalb/'^X^ below which af^ = af = = and 
— /^2* — ^/ih- Above this threshold, we find 



^2 = ± A/^(e-eth) 



= = ±y^(e-eth), 

/?r=/5r = ^, (43) 

where all the aj must have the same sign. What is unusual about this system in 
comparison with the normal optical parametric oscillator is that it is stable at threshold, 
with the critical point for the below threshold solutions being at a pump amplitude of 
Cc = lalb/V^X^ which means that ec = V^eth- For this pumping strength, the actual 
above threshold solutions are stable, so that this system, at least for the parameter 
regimes we consider here, is always stable. With these steady-state solutions, we have 
all the information we need to calculate the output spectral correlations of interest. 
Below the oscillation threshold we can find relatively simple expressions for the 
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Figure 10. Tripartite entanglement criteria for the system of section 14.21 for 
7a = 76 = 1, X = 10"^ and e = 2eth. 



output spectral correlations, 



2^2 



Si^{uj) 



5- 



(44) 



whereas above threshold the analytical expressions become more complicated. 

We show results for the tripartite criteria in figure IHl immediately below the 
oscillation threshold, and in figure ^1 for e = 2eth- We find that the tripartite 
entanglement is strongest at threshold, with = >S'23 and 5*13 showing a lesser violation 
of the criterion. We also found that, well above threshold, the individual quadratures 
Xi and X3 become squeezed, although X2 remains above the shot-noise level. 



5. Conclusions 

We have examined two different interaction schemes in terms of their potential for 
creating continuous variable tripartite entanglement. One of these is based on two 
cascaded nonlinearities and the other on two concurrent nonlinearities. Although the 
scheme based on cascaded linearities has been analysed previously, we have extended 
these analyses to include depletion of the pumping fields. Both these schemes were 
shown to exhibit approximately the same degree of tripartite entanglement, and are 
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therefore suitable candidates for practical applications. The concurrent scheme has the 
possible advantage that two of the outputs are of equal intensity and is symmetric in two 
of the correlations, whereas the cascaded scheme produces three different intensities and 
all three correlations give different values. Which one is preferable in a given situation 
may come down to the preferences and experiences of the experimenters wishing to use 
them. 
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